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Inevitable Irreversibility Generated by the Glass Transition 
of the Binary Lattice Gas Model 
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We numerically investigate the thermodynamic properties of the glass state. As the object of 
our study, we employ a binary lattice gas model. Through Monte Carlo simulations, we find that 
this model actually experiences a glass transition. We introduce a potential into the model that 
represents a piston with which we compress the glass. By measuring the work performed in this 
process, we find that irreversible works exist at the glass state even in the quasistatic limit. This 
implies that yield stress is created by the glass transition. 

PACS numbers: 61.43.Fs 65.60.+a 
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Theoretically describing structural glass is a long- 
standing problem in statistical physics. As the tempera- 
ture of a supercooled liquid is reduced, its molecular dy- 
namics become slower, and eventually such motion ceases 
at a temperature T g , which is called the glass transition 
temperature |l], g, f| |j. Below this temperature, the 
glass is usually considered to be a solid in the sense that 
its molecules are trapped within a cage formed by their 
neighbors from which they cannot escape, at least within 
the experimental time scale. However, a glass is never a 
solid in the true sense, but out of equilibrium. 

Significant effort has been devoted to understanding 
the nature of glasses, with theoretical studies employing 
mode coupling theory || , the replica model |6[ fjj , and the 
trap models |g, However, in spite of these attempts, a 
unified picture and satisfactory understanding have not 
yet been obtained. The main source of the difficulty in 
obtaining a theoretical description is that we cannot use 
the universal measure (i.e. the canonical measure) in 
treating the glass state. Since the glass state is not an 
equilibrium state, the Boltzmann-Gibbs ensemble does 
not apply. Without the framework of statistical mechan- 
ics, we are powerless to theoretically treat a large number 
of degrees of freedom. To allow for the treatment of glass 
forming materials, it is thus necessary to construct an ex- 
tended form of statistical mechanics. With this goal, as a 
first step, it is wise to first consider the thermodynamics 
of glasses. 

Since a glass represents a nonequilibrium state, it is 
not definite whether thermodynamic structure exists in 
this state. For this reason we first reconsider the thermo- 
dynamic concept of a glass state. After confirming the 
important macroscopic nature, we hope to connect the 
macroscopic picture to the microscopic picture. In this 
letter, we report the numerical investigation of the ther- 
modynamic properties of the glass state in a lattice gas 
model. 

Lattice gas models are very simple and easy to treat 
computationally. For this reason, if a lattice gas model 
possesses a glass state, it would be convenient for the 
theoretical study of this state. Several lattice gas models 
have been proposed as models of the glass transition. The 



most commonly used models are the lattice gas models 
with kinetic constraints proposed by Jackel, Palmer [Q, 
Kob and Andersen and Kurchan et al., among oth- 
ers. However, because these models possess no energy, 
then are not suited to our thermodynamic study. Despite 
this fact, we are interested in lattice-gas models because 
Fuchizaki and Kawasaki assert that such a model with 
a repulsive interaction displays glassy dynamics Jl2| p"3| . 
Since a lattice gas model of this type does possess an en- 
ergy, it would seem to be a possible model with which to 
carry out our study. However, in our simulations, we have 
not found a glass transition in this model. (We give de- 
tailed discussion regarding this point below.) Therefore, 
we must make a new lattice model with energetic con- 
straints. In this letter, we extend the lattice gas model 
to a two component model. Using Monte Carlo simula- 
tions, we find that our extended lattice gas model indeed 
exhibit a glass transition. Finally, we carry out thermo- 
dynamic operations within the model and investigate its 
macroscopic properties so revealed. 

We first introduce our extended lattice gas model. We 
consider a system defined on a 2-dimensional lattice of 
size L x L. We assume that there are two types of parti- 
cle, A and B, on the lattice, and that all particles inter- 
act repulsively. The strengths of the interactions differ 
for different pairs, A and A, B and B, and A and B. 
The interactions are restricted to nearest neighbors, and 
there is no more than one particle at each site. These 
conditions are satisfied by the Hamiltonian 
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where nf — 1 and nf = if site i is occupied by a particle 
of type A, while nf = and nf = 1 if it is occupied by 
a particle of type B, with i = 1,2,- • • L x L labeling the 
L x L lattice points. Here, (ij) denotes nearest neighbor- 
pairs. The quantities J AA , J BB , J AB are the coupling 
constants between A and A, B and B, and A and B. We 
call this model the "binary lattice gas model" . The model 
of Fuchizaki and Kawasaki is single component version of 
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FIG. 1: Monte Carlo simulation results for the density- 
density correlation function <f)(t + t w ,t w ) of the binary lat- 
tice gas. Correlations are plotted for several waiting times 
t w : t w = 0, 10 3 , 10 4 , 10 5 , 10 6 . (a) The parameter values are 
here (N A ,N B 1 T) = (65,65,1000). The unique curve in 
the long waiting time limit has a finite Debye- Waller fac- 
tor fn defined by fo = limt_> 00 lim tl „_>oo <$>■ The nature of 
the correlation dynamics in this case indicates that the sys- 
tem is in the glass state, (b) The parameter values are here 
(N A , N B ,T) = (55,55,1000). In this case, the unique curve 
in the long waiting time limit has a vanishing Debye- Waller 
factor. The nature of the correlation dynamics in this case 
indicates that the system is in the gas state. 



when the particle concentration is sufficiently high. In 
Fig. 1 (a) we show the results for the density-density cor- 
relation functions with several waiting times for the case 
T = 1000 and N A ,N B = 65. We see that the correla- 
tion functions exhibit an aging dynamics, that is, the 
correlation depends explicitly on the observation time 
and on the time at which the perturbation is applied. 
We note that when the waiting time is sufficiently long, 
the correlation function converges to a universal form, 
which has a finite Debye- Waller factor. Thus we con- 
clude that the system is in the glass state in this case 
with N A ,N B = 65. In Fig.l (b), we show the density- 
density correlation functions for the case T = 1000 and 
N A ,N B = 55. We see that again in this case, the cor- 
relation function converges to a universal form, but here 
it has a vanishing Debye- Waller factor, so that the sys- 
tem is concluded to be in the gas state in this case with 
N A ,N B = 55. We can hence conclude that the tran- 
sition between the gas state and the glass state occurs 
at a concentration in the range C = 0.42 — 0.51. The 
precise determination of the transition point is difficult, 
because the relaxation becomes very slow near the transi- 
tion point. However, fortunately a precise determination 
is not necessary for the thermodynamic study we wish 
to carry out. For this purpose, it is only important that 
our binary lattice gas model has a concentration glass 
transition. (We have found that this model has another 
type of transition, a glass transition as a function of tem- 
perature. However, this transition is not involved in the 
present work.) 

Now we investigate the thermodynamics of the model. 
For this purpose, it is convenient to slightly modify the 
model described above. First, we change the boundary 
conditions. Among the four boundaries i x = 0,i y = 
0, i x = L + 1 and i y = L + 1, we consider the pair of 
boundaries i y = and i y = L + 1 to be still periodic, 
while we consider the other pair i x = and i x = L+ 1 to 
be perfectly impermeable. Next, we introduce a potential 
which represents a movable piston into the model. The 
piston potential, denoted by V a (i), is added to the origi- 
nal, "bare" Hamiltonian Ho, so that the full Hamiltonian 
Ji is given by 



this model. We stipulate that the dynamics of the system 
are conservative by enforcing the conditions 

Y^nf=N A ^nf =N B . (2) 

i i 

To enforce these conditions, we adopted the Kawasaki 
(spin-exchange) dynamics for the numerical simula- 
tions jlB). Thus, the system is parametrized by the set of 
quantities J AA , J BB , J AB , N A ,N B ,T and L, where T is 
the temperature of the system. In our simulations, nega- 
tive coupling constants were used in order to preclude 
the gas-liquid transition. We set ( J AA , J BB , J AB ) = 
(— 1, —4, —2) and L = 15, and employed periodic bound- 
ary conditions. The results of our Monte Carlo simu- 
lations indicate that this model possesses a glass phase 



H = H + J2 v a(*)( n ?+ n ?)- ( 3 ) 

i 

We represent the piston as a quadric potential: 

VS) = [^-f (*.>■) (4) 

[0 (i x < a). 

By changing the value of a, we carry out the operation 
of compressing or expanding the system, and we refer to 
its value as the "piston position". In our numerical sim- 
ulations, a was changed as follows. The piston position 
is initially L, and at this time, the system is in the gas 
state. Then we compress the system gradually until a 
is Lq. It is then gradually returned to L. We used the 
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following simple form of a(t): 



o(t) = 



-vt + L 
v{t-T) + L 



(0 < t < T) 
(T <t< 2T), 



(5) 



where v is a parameter that represents the speed of the 
operation and T is defined by T = (L — L )/v. When 
we make Lq sufficiently small, the system experiences 
the glass transition during the compression. With the 
transition between the gas and glass states caused in this 
manner, we investigate the thermodynamic properties of 
the glass state. The quantity that we measure in the 
simulation is the work performed on the system, defined 
by 



! )da. 



(6) 



We are particularly interested in the quasistatic limit of 
this function, and we write this limiting form as W qs . 
In equilibrium, W qs is equal to the Hclmholtz free en- 
ergy difference between the initial and final configura- 
tions. However, because a glass state is non-equilibrium 
in nature, it is not certain if W qs also can be used to 
construct the free energy in this case. (Here, we note 
that we use the word "quasistatic" to mean simply that 
the process is carried out sufficiently slowly compared to 
the /3-relaxation time scale.) The main question we wish 
to answer in this letter is whether W qs is a useful ther- 
modynamic function in the description of the glass state. 
For this purpose, we investigate the v dependence of the 
work W[a] in a complete cycle of compression and ex- 
pansion defined by (g). If lim„->o W[a] = in this cycle, 
Wq S constitutes a thermodynamic function of the glass 
state, as it is in equilibrium thermodynamics. 

In this simulation we set the parameters as L = 15, 
T = 1000 and N A ,N B = 30. Lq is taken sufficiently 
small as the system becomes a dense packing state at 
a = Lq. With the parameter values given above, the 
glass transition occurs at approximately a — 9. In Fig. 2 
(a) we display simulation results for W[a] as a function 
of v. For comparison, in the same figure, we also plot the 
results obtained using a Fuchizaki-Kawasaki type system 
with N A = 60, N B = 0. As we see, in this case the 
system always remains in the gas state. We can clearly 
see the effect of the glass transition from this compari- 
son. For each case, we see that the work converges to a 
finite value in the quasistatic limit, but we see that the 
irreversible work performed in the two component case 
is much larger than that in the single component case. 
In the single component case, the amount of irreversible 
work ins very small, and it appears to be due to a lat- 
tice effect. This is demonstrated by the fact that the 
amount of this work does not change when we set the 
values of coupling constants to zero. Thus we presume a 
slight amount of the work of the single component case 
to be almost zero. Contrastingly, in the two component 



case, irreversible work is essential to the glass transition. 
This implies that there is no thermodynamic function on 
characterizing the glass state. 

The effect of the glass transition is seen more explicitly 
when we consider the dependence of the quasistatic irre- 
versible work on the system size. In Fig. 2 (b) we display 
the simulation results for W[a] as a function of v for two 
systems with four times more lattice points, but the same 
density: L — 30 and (N A ,N B ) = (120, 120) for the two 
component case, and L — 30 and (N A ,N B ) = (240,0) 
for the single component case. We find larger difference 
here between the value of W qs for the one component 
and two component systems than in the L = 15 case. 
In Fig. 3 we display the system size dependence of W qs 
for the two component system, with equal densities in 
each case. We find that the the irreversible work exhibits 
an anomalous and nonextensive system size dependence: 
Wqs ~ (L x L) 2 . The reason for this dependence is not 
yet clear. 

From these results, we conclude that the glass state 
does not conform to conventional thermodynamics laws. 
We cannot define a thermodynamic function describing 
the glass state, and we need a theory which deals with 
the inevitable irreversibility directly. We point out here 
that approximately a decade ago, Langer and Senthna 
reported that the glass transition exhibits a hysteresis of 
the entropy |ll| . Our results may be considered the free 
energy version of their results. 

We have arrived at an interesting conclusion. From 
the operational point of view, the glass transition can be 
described as a generation of irreversibility, even in the 
quasistatic limit. For the description of a glass, it is thus 
not possible to use a simple extension of conventional 
thermodynamics. It is interesting that the macroscopic 
nature of our lattice gas may be represented by the transi- 
tion between a Newtonian fluid and a Bingham fluid p7| , 
because the inevitable generation of irreversibility in our 
system can be regarded as a generation of yield stress. In 
fact, we can determine the yield stress from Fig. 2. How- 
ever, at present we do not know of a clear relation be- 
tween a glass and a Bingham fluid. The most important 
properties of glasses are disorder and metastability, and 
it is possible that these properties correspond to prop- 
erties of Bingham fluids, but this correspondence is yet 
to be found. In any case, it may be productive to recog- 
nize that the glass transition may be one of the route to 
Bigham nature. 
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FIG. 2: The irreversible work as a function of the speed of 
the compression-expansion operation for both lattice gases. 
The speed of operation v is measured in units of the inverse 
of the time needed to change the system length by the inter- 
site distance, (a) The parameter values are here L = 15, 
T = 1000, and the particle numbers are (N A ,N B ) = (60,0) 
and (N A ,N B ) = (30,30). For each case, the amount of 
work converges to a finite value in the quasistatic limit. The 
amount of quasistatic irreversible work in the two compo- 
nent case is much larger than that in the single compo- 
nent case. (b) Here the parameter values are L = 30, 
T = 1000, and the particle numbers are (N A ,N B ) = (240,0) 
and (N A ,N B ) = (120, 120). We see a more distinct difference 
between the two component and single component cases here 
than in (a). 
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FIG. 3: Log- log plot of the system size dependence of the 
quasistatic irreversible work in the two component case. The 
density is the same in each case. System size is normalized as 
15 x 15 is equal to 1. The dashed line represents the fitting 
function W ~ (system size) 2 . 



